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Abstract 

We argue that the T nuclear modification factor, /?J^^, measured at RHIC can only 
be reproduced once gluon EMC effect and fractional energy loss are taken into account. 
At backward rapidities, the visible suppression of hints at the presence of a gluon 
EMC effect, analogous to the quark EMC effect - but likely stronger. At forward and 
mid rapidities, the data can only be accounted for by fractional energy loss, recently 
revived in the literature. Our conclusions do not depend on the detail of the nuclear parton 
distributions. Such a fractional energy loss also provides an alternative explanation to gluon 
saturation for the strong forward Jj^ suppression in dka collisions observed by PHENIX. 



1 Introduction 

Conventional nuclear modifications of gluon distribution in heavy ions - known as shadowing 
and anti- shadowing - as well as the possible break up of the bb pair along its way off the 
nucleus are shown to have a limited impact on the T production in JAu collisions at RHIC at 
^fsNN = 200 GeV. Without additional effects, one fails to reproduce the nuclear modification of 
the yield measured by the STAR and PHENIX experiments [1,2]. 

This discrepancy motivated us to study the impact of two additional nuclear effects. The first 
is a suppression of gluon distribution in nuclei for intermediate Bjorken-;c, 0.35 < Xb < 0.7, 
analogous to the quark EMC effect [3], unobserved until now. The second is di. fractional energy 
loss [4] proportional to the projectile parton energy and caused by medium-induced radiations 
associated to the quarkonium hadroproduction. This radiative energy loss arises when radiation 
off the incoming parton and outgoing coloured object is nearly coherent, which happens when 
the flow of color charge undergoes a small angle scattering in the nucleus rest frame. Such 
fractional energy loss is thus still relevant at RHIC energies, contrary to the usual radiative 
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energy loss [5] where a hard parton is suddenly produced (or undergoes large angle scattering) 
in a static medium. 

We have compared our results -from our Monte-Carlo framework JIN [6]- to dAu Y data 
from RHIC experiments [1,2] and found a good agreement only when including both additional 
effects. We also found a good agreement with dAu J/i^ data when the effect of energy loss 
is combined with mild gluon shadowing and a small value for the effective absorption cross 
section. In particular, the energy loss provides an alternative explanation to gluon saturation 
-or strong shadowing- for the strong forward J/ilf suppression in dAu collisions observed by 
PHENIX. Whereas, in the T case, energy loss seems to be required since no other effect can 
explain the forward data, energy loss is only one of the possible effects responsible for forward 
7/(/f suppression. The situation is indeed less clear for 7/i/r since shadowing and/or saturation 
as well as the break-up probability can also have a significant impact [6]. 

2 Propagation in cold nuclear matter: T\s J/iJ/ 

The first effect to be discussed is the probability for the heavy-quark pair to survive the propa- 
gation through the nuclear medium, usually parametrised by an effective cross section cTes. It is 
sometimes referred to as the nuclear absorption or break-up probability. 

A priori, the smaller size of the bb pair when compared to cc pair implies that bb states 
should suffer less break-up than cc. Yet, the ratio of their size depends on the evolution stage 
of the heavy-quark pair. At the production time, this ratio is expected to be mh/nic, hence a 
size 3 times smaller for bb than for cc. When they are fully formed, it is rather ^4t^ x — 
as expected from their Bohr radii, hence a size 2 times smaller. The relevant timescale for the 
pair evolution is the formation time. According to the uncertainty principle, it is related to the 
time needed - in their rest frame - to distinguish the energy levels of the 15 and 25 states [7]: 
tf = =2x10 GeV / 10.5 GeV^ = 0.4 fm for the T. 

For our purpose, tf has to be considered in the rest frame of the target nucleus, i.e. the Au 
beam at RHIC. The relevant y factor is then obtained from the rapidity of the pair corrected by 
the Au beam rapidity y = cosh(j - y^^^^) where J^^",,, = -5.36 for RHIC. The formation time 
for the different rapidities reached by RHIC experiments are given in Table. 1 . tf is significantly 
larger than the Au radius - except in the most backward region. This implies that the bb is 
nearly always still in a pre -resonant state when traversing the nuclear matter. 

Table 1: Y boost and formation time in the gold rest frame as a function of its rapidity at 
= 200 GeV. 
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For the forward and mid rapidity regions, this has two consequences : first, the break-up 
probability is expected to be small, of the order of a tenth of that of J/i/r {{mdnibf ~ 0.1), 
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following the early-time scaling mh/nic; and second, it ought to be the same for the TCli"), 
T(25 ) and Y(35 ) states, since these states cannot be distinguished at the time they traverse the 
nucleus. At the LHC, t/s are even larger; excited Y state suppression in PbPb collisions seen 
by CMS [8] can thus only be explained by hot nuclear effects. 

The backward region at RHIC requires a closer look. Indeed, for the same tf, of the order of 
15 fm, the E866 experiment at Fermilab [9] observed a different suppression of J/i^ and if/(2S) 
produced with Feynman Xf up to 0.2 at yfsNN = 38.8 GeV. We might thus expect different 
absorption cross-section for the T(1S), Y(25) and Y(35) states for y < -1. However, the E772 
experiment at Fermilab [10] has measured the Y(15) and Y(25 -I- 35) separately at = 38.8 
GeV down to negative Xf - with even smaller tf - and it observed a similar suppression for the 
IS and the (2S +3S) states The only explanation for such a result is that the absorption of the 
bottomonium is actually very small, preventing us to see any difference of absorption between 
the 3 states. In the following, we will consider a range of cTes from to 1 mb even though 1 mb 
has to be seen as a conservative upper bound. 

3 Gluon momentum distribution in nucleus 

At high energy (small xb), the nuclear Parton Distribution Functions (nPDF) differ from those 
of free nucleons due to non-linear QCD effects. Nucleons shadow [11, 12] each other and the 
nPDFs are expected to be lower than for free nucleons. At 0.01 < xg < 0.3, some experimental 
data hint [13] at an excess of partons with regards to unbound nucleons, referred to as anti- 
shadowing. For 0.35 < xb < 0.7, the distribution is depleted again. This suppression is known 
as the EMC effect. 

These nuclear modifications are usually expressed in terms of the ratios Rf of the nPDF 
of a nucleon bound in a nucleus A to the free nucleon PDF. The numerical parametrisation of 
Rfi^B, Q^) is usually given for all parton flavours. Here, we limit our study to gluons since, at 
RHIC, Y is essentially produced through gluon fusion [14, 15]. To best explore the possible 
impact of we have considered 3 sets: EKS98 [16], EPS08 [17] and nDSg [18] at LO. 
Recently, a new set with fit uncertainties, EPS09 [19], has been made available. Yet, in the case 
of gluons, nDS and EPS08 match -except for xb > 0.3- the extreme values of EPS09. Besides 
EKS98 is very close to EPS09 central values. We thus consider more illustrative to use EKS98, 
EPS08 and nDSg. The spatial dependence of PDF nuclear modification has been included with 
a modification proportional to the local density [20] . Following the common practice, we label 
X\ {X2) the gluon momentum fraction in the proton/deuteron (nucleus). 

Thirty years ago, the European Muon Collaboration (EMC) [3] observed a depletion of the 
quark densities in nucleons bound in nuclei, when compared to the ones of free nucleons, in 
the range of 0.35 < xb < 0.7. To date, there is no consensus [21] about the origin of this 
suppression, referred to as the EMC effect. It is still the subject of vivid activities. It has been 
attributed to local nuclear density effects, to properties of the bulk nuclear system, and recently 
[22] to Short Range Correlations (SRC) in the nucleus. Up to now, this effect has not been 
confirmed for gluons, even if it is allowed in some of the shadowing fits. 

While gluon shadowing in the existing nPDF constrained fits, especially in EPS 08 & 09, 
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Figure 1: The 3 gluon nPDF sets we have used to single out the EMC effect and the correspond- 
ing EPS 09 uncertainty for Au. 

is the subject of intense on-going debates, the gluon EMC suppression is usually overlooked. 
Indeed, very little is known about gluons in this region and few data constrain their distribution 
at xb larger than 0.3. The amount of the EMC suppression is actually pretty much unknown [19], 
except for a loose constrain set by momentum conservation. We also note that another effect 
[23] arising from momentum conservation within the nucleons could be at play at larger xb, 
Xb > 0.7, i.e. at larger y or smaller ^/snn, and hence it is not applicable here. 

To single out a possible impact of the EMC gluon suppression on T, we have used three 
of the EPS09 LO sets: one with a quark-like EMC gluon suppression, and the two limiting 
curves in the region 0.35 < xb < 0.7 (Fig. 1). As the data-theory comparison will show, its 
magnitude indeed seems stronger than what has been previously supposed in the existing nPDF 
sets (which assumed a quark-like EMC gluon suppression) and disagrees with a Fermi-motion 
enhancement down to Xb - 0.5 as expected in [24]. 

In practice, to account for the nuclear effects on Y production in nucleus collisions, we use 
our Monte-Carlo framework JIN [6], based on the probabilistic Glauber model, used to describe 
7/i/^ production at RHIC [25, 26]. It allows to consider improved kinematics corresponding to a 
2 ^ 2 (g + g ^ bb + g) partonic process for the T production (as in the Colour- Singlet Model 
(CSM) at LO [27]). In earlier studies of nuclear matter effects on Y production [28], the bb pair 
was assumed to be produced by a 2 ^ 1 partonic process (^ -i- ^ ^ bb). 



4 Energy loss 

An energetic parton travelling in a large nuclear medium undergoes multiple elastic scatterings, 
which induce gluon radiation, hence a radiative energy loss. Effectively, since the parton loses 
energy, the quarkonium hadroproduction process probes the incoming proton PDF at a higher 
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x[ = Xi+ Axi, leading to nuclear suppression since the gluon PDF g(x) atx = x\ is smaller than 
the one at x = Xi, g(x'^ < g(xi). 

Twenty years ago, parton energy loss in nuclear matter was suggested to be the dominant 
effect responsible for both Drell-Yan and J/i/r suppression [29] at forward Xf. It was assumed 
to scale with the incoming parton energy, AE oc E, while being a medium-induced energy loss 
of an independent colour charge taking place in the medium and not in the vacuum. However, 
general quantum mechanical arguments based on the uncertainty principle seem to lead to the 
bound AE < L{k\) [30], where L is the size of the medium and {k\) is the average transverse 
momentum squared of the radiated gluon. Such medium-induced energy loss therefore cannot 
scale with the parton energy E when E ^ oo. 

Recently [4], it has nevertheless been pointed out that, in the case of small angle quarkonium 
production, the medium-induced spectrum also arises from large formation times tf » L. The 
uncertainty principle does not constrain anymore the radiation and the medium-induced "energy 
loss" scales as the energy E in apparent contradiction with [30]. The contradiction is solved by 
noting that the "loss" derived in [4] is not the loss of an independent parton suddenly produced, 
but the amount of gluon radiation associated to the quarkonium production process. This applies 
for octet-like mechanisms where the heavy-quark pair is produced at short distances in a colour- 
octet state - the latter can thus radiate. This is so in the Colour-Evaporation Model (CEM) and 
the Color-Octet Mechanism (COM) (see [14, 15]). For singlet-like mechanisms - such as the 
CSM - the energy loss will occur only when the heavy-quark pair is produced at short distances 
(in a colour-singlet state) in conjunction with at least one hard gluon (or quark). This is the case 
in our CSM LO 2^2(^-i-^^T-l-^) process valid at low Pt [27], but also at higher-orders 
in as with 2 or 3 coloured partons in the final state [31, 32]. 

One then finds [4] that for forward angle T production in pA collisions, the fraction of 

medium-induced radiated energy is given by AEfE = Ax^lx^ - NcU^ ^Ja{p^)/Mt, where 

A(/7y) is the broadening of the radiated gluon from the proton and Mj is the transverse mass of 
the final-state coloured object. We assume the gluon broadening to be equal to the Y broadening, 
A{pj) = {pj){A) - (p|)(^H), proportional to the length of nuclear matter seen by the incoming 
parton. This broadening in p^ is consistent with a dependence on A^^^, as expected in multiple 
scattering models. Taking the E772 value measured with W target [33, 10], A{p^) = 0.410 
GeV2, one can thus write A{p^^) = 0.072 GeV^ A^'\ 

In octet-like mechanisms, Mj is the mass of the quarkonium bound state, here Mj ~ 10 
GeV. For singlet-like production, Mr is rather the mean pr of the final-state gluon. We have 
taken Mj = 3.0 GeV here. For Nc = 3 and as = 0.2, we have Ax"'^^'"' lx\ ~ 13% in the singlet 
case, and about three times lower, Ax^['^^"'' I x\ ~ 4%, in the octet case. As we shall see, this 
makes a visible difference in the yield suppression. At yfsNN = 200 GeV, the average xi of the 
projectile gluon in the forward rapidity region, 1.2 < y < 2.2, is = 0.28. In the mid 

rapidity region, (x™'''^) = 0.05. The corresponding suppression is just the ratio between the 
PDF of the gluon evaluated at its original x\, x\ = xi + Ax\ before the loss and the PDF at the 
resulting x^: RhAxi, Q^) = g(x[, Q^)/g(xu Q^). 

Before evaluating the corresponding yield suppression, it is important to note that, in the 
CSM, Xbi production at LO does not produce a recoiling gluon and is not affected by such an 
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effect. We evaluate it to be 10 % of the Y yield. In the other cases, the feed downs are expected 
to be suppressed as much as the direct yield in the corresponding model. Taking modem global- 
fit gluon PDFs, one finds that this fractional energy loss gives R]"^^'"^' = 65% 70% and ^^^'f = 
85% 90% at forward rapidity. At mid rapidity, there is still a suppression due to this energy 
loss effect, which induces Rf^f^''" ~ 80% and R^ffJ ~ 90%. We have recomputed these numbers 
using different modern LO PDF sets and have found that our result lies well in this range. We 
consider the overall uncertainty on 1 - to be of the order 10-20 %, i.e. clearly smaller than 
that of the data. 

We have also studied the impact of the energy loss effect on J/tff production. Due to the 
smallness of Mr compared to the T case, we have Ax-f"^'"/xi ~ 20% and Axf*7xi ~ 10%. 
Both losses are stronger than for the Y. Yet, at ^/snn = 200 GeV, {x^^'"^'^^^) ~ 0.1 and 
(^^midji^^ ~ 0.02, i.e. smaller than for Y. Hence, one has i?f"f' ~ 70% and R';^' ~ 80% 
for forward J jifj. At mid rapidity, the induced suppression is R""f^''" ~ 80% and R"^'^''^' ~ 90%, 
similarly to the Y case. These should however be combined with the suppression from the 
shadowing and nuclear breakup which -contrary to the Y case- are not negligible. 

At backward rapidities, the typical formation times are of the order of the nucleus size and 
the conditions of the applicability of the derivation of AE in [4] are not satisfied. In addition, 
due to the behaviour of the gluon PDFs at small xb, one would find a smaller effect than at 
J > 0, if one were tempted to apply the same formula without much concern. 



5 Results for dAu collisions 

Experimentally, Y and J/iJ/ suppressions are studied by measuring a nuclear modification factor 
RdAu, the ratio of the yield in dAu collisions to the yield in pp collisions at the same energy 
times the average number of binary inelastic nucleon-nucleon collisions Ncoii in a dAu collision: 
^rfAu = {dNdAu)/{{Ncon)dNpp). Any nuclear effect affecting Y and J /iff production leads to 

From STAR and PHENIX [1,2] only the rapidity dependence is known. For now, the 3 
Y resonances are not resolved but are measured together. Since the nuclear absorption has to 
be small and since the nPDF effects are very likely similar for these 3 states, we can safely 
consider them on the same footage. However, it is worth noting that a future measurement of 
RdAu focusing only on YCl^) would be very precious to confirm this assumption. 

Fig. 2 (a) shows the uncertainties on R^au vs y due to the lack of knowledge on gluon nPDF 
and due to a variation of cr^g between and 1 mb. Gluon shadowing (at forward y) and anti- 
shadowing (at mid y) are not strong enough to match the RHIC data. The nuclear absorption 
also shows a very mild effect, which is insufficient to reproduce the data, especially since 1 mb 
is likely exaggerated. At backward y, a gluon EMC effect stronger than that of quarks, such as 
the one provided by the EPS09 lower bound, gives a very convincing account of the backward 
data (see Fig. 2 (b)). This modification does not alter the results at mid and forward rapidities 
and is perfectly legitimate given the current knowledge of the gluon nPDF in this region. It is 
nevertheless clear that experimental uncertainties are large and the observation of such a strong 
EMC gluonic effect would only be confirmed once we have data with reduced errors. The data 
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■ Production-kinematics ^ Singlet-like 
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Figure 2: (b) Theoretical uncertainties on R^^^^ due to the nPDF (coloured lines) and oss (purple 
band).; (c) Effect of an increase of the gluon EMC effect and of the fractional energy loss for 
singlet-like (red band) and octet-like (green band) production at small angles {y > 0) on i^J^^. 
(d) i?^^ from [26, 34] (top blue band) for (Tabs = 2 mb with the fractional energy loss for 
singlet-like (red band) and octet-like (green band) production at small angles. Data for T are 
from [1,2] and for 7/i/^ from [36]. 

however already visibly disfavour the case with no gluon EMC effect (orange solid bars). 

We now discuss the effect of the fractional energy loss. Fig. 2 (b) clearly shows its impact in 
the mid and forward regions. The green bands indicate the suppression expected for octet-like 
productions, and the red bands for a singlet-like production. The second case - with a stronger 
suppression - better matches the data which cannot otherwise be accounted for. Going further 
in the interpretation, we could say that the significant suppression of T in the forward region 
favours a production mechanism where the heavy-quark pair is produced at short distances in a 
colour-singlet state in association with one or more gluons. 

Such energy loss also impacts 7/i/r production, as shown Fig. 2 (c), where our previous 
results with EKS98 [26, 34], for an ad hoc (Tabs = 2 mb, are combined with the additional 
suppression induced by the energy loss. Our results are qualitatively confirmed by a recent 
analysis [35] where energy loss is also shown to reproduce the most forward fixed target data 
at lower energies. Agreement with PHENIX data [36] is obtained especially if one assumes 
(Tabs to be negligible as done in [35]. In particular, the suppression due to the energy loss has a 
steeper y dependence than that of shadowing. 

6 Conclusion 

We have investigated T suppression in dAu collisions at ^/s^n = 200 GeV. The rapidities cov- 
ered by RHIC experiments allow for a unique study of cold nuclear matter effects and revealed 
unexpected features presented here. Backward rapidities correspond to the largest xj, above 0.2, 
where one expects an EMC suppression - at least for quark PDF. For mid and slightly backward 
rapidities, one expects anti- shadowing, i.e. an excess of partons inducing an excess of Y. Fi- 
nally, the forward domain - where X2 becomes small - should be subject to parton shadowing. 
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giving a reduction of the yield. 

Our findings are as follows: in the most backward region, studied by PHENIX, the sup- 
pression of the T yield in the dAu data may be a first hint that the EMC suppression of the 
gluon distribution is stronger than the quark one. It is probably the very first observation of 
such an eff"ect, whose quantitative understanding may in the future provide us with fundamental 
information on the internal dynamics of heavy nuclei. 

Due to the large scale set in by the Y mass, shadowing is found too feeble to reproduce data 
at y > 0. Moreover, the data at 3^ ^ does not exhibit any excess which would pin down anti- 
shadowing. This motivated our study of the recently revived medium-induced fractional energy 
loss, this mechanism which would also explain the anomalous suppression of the forward J/il/ 
at E866 and NAB in pA [4]. 

We have shown that -at variance with the J/i/r case [26]- gluon shadowing has, for any 
nPDF fits, a small effect at the X2 and of forward-T production at RHIC. An X2 scaling of 
the nuclear effects is thus not expected. Provided that (fractional) energy loss is the dominant 
effect at work in the forward region, an xi scaling is expected, which exhibits an xp scaling in 
the forward region where Xi - Xf,'m analogy with [37, 10]. 

For bottomonia, whose survival probability to escape the Au nucleus should be large and 
whose shadowing suppression is found small, we argued that a new mechanism such as this 
fractional energy loss seems to be required. In the J/tff case, it is rather allowed. The same 
mechanism may still apply for y ^ - still in the small angles domain at RHIC - and also 
elegantly explains the trend of the T STAR data. Such an effect provides an alternative expla- 
nation for the strong forward-7/i/^ suppression in dAu as seen by PHENIX [36]. Our results are 
confirmed in [35]. Quantitatively, Y data are better explained if the bb pair is produced at short 
distances in a colour-singlet state with a hard gluon as in [27], as opposed to a production in a 
colour-octet state. The picture is less clear for J/ifr. 
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